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! We present models for the HI 21 cm emission from circumstellar disks and 

use them to convert observed upper limits on the Hi 21 cm flux to limits on 
.£h ! the total disk gas mass. The upper limits we use come from earlier Australia 

Telescope Compact Array observations of the debris disk around (3 Pictoris as well 
as fresh Giant Meterwave Radio Telescope observations of HD 135344, LkCall5 
and HD 163296. Our observations and models span a range of disk types, from 
young proto-planetary disks to old debris disks. The models self-consistently 
calculate the gas chemistry (H/H2 balance) and the thermal structure of UV 
irradiated disks. Atomic hydrogen production is dominated by UV irradiation 
in transition phase objects as well as debris disks, but for very young disks, HI 
production by stellar X-rays (which we do not account for) is important. We 
use a simple radiative transfer approach to convert the model disk parameters 
into predicted Hi 21 cm line maps and spectral profiles. This allows a direct 
comparison of the observations to the model. We find that the Hi traces the disk 
surface layers, and that the H I emission, if detected, could be used to study the 
effects of irradiation and evaporation, in addition to the kinematics of the disk. 
Our models cover massive protoplanetary disks, transition phase disks and dusty 
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debris disks. In massive protoplanetary disks, UV produced H I constitutes less 
than 0.5% of the total disk mass, while X-rays clearly dominate the chemistry 
and thus the Hi production. For the two such disks that we have observed, 
viz. those around LkCa 15 and HD 163296, the predicted 21 cm flux is below the 
current detection limit. On the other hand, transition phase disks at distances of 
100 pc have predicted 21 cm fluxes that are close to the detection limit. Finally, 
in debris disks, hydrogen is mainly molecular since the high dust-to-gas mass 
ratio leads to warmer disks, thus increasing the formation rate of H 2 . This, in 
conjunction with the small total gas mass, makes the predicted flux to fall below 
our detection limit. However, future telescopes like the SKA should be able to 
image the Hi 21 cm emission from nearby transition phase disks, while a radio 
telescope with only ~ 10% the area of the SKA should be able to detect the 
emission from such disks. Such 21 cm line observations would probe both disk 
evaporation as well as disk kinematics. 

Subject headings: accretion, accretion disks - circumstellar matter - stars: for- 
mation, pre-main-sequence - infrared: stars 

1. Introduction 

Protoplanetary disk evolution and planet formation are closely intertwined. To under- 
stand planet formation, we first need to develop a comprehensive picture of the physical and 
chemical conditions in protoplanetary disks. While dust has been observed over the entire 
range of disk evolution — from the very young massive disks out to the old debris disks that 
contain only a few lunar masses of small dust grains — , gas is much more difficult to detect. 
Hence our picture of the chemical evolution of the gas phase is very incomplete and mainly 
restricted to a few prototypes. This is unfortunate, because gas plays an important role in 
the disk physics. It determines the hydrostatic equilibrium in the early phases, it affects the 
dust grain growth and dynamics as well as planet migration, and eventually its dispersal 
ends the planet formation process. 

Hydrogen is the most abundant element in molecular clouds and in the protoplanetary 
disks that are a byproduct of the star formation process. Most of the disk mass is in the form 
of either atomic or molecular hydrogen. Depending on the number of photons capable of 
ionizing hydrogen (i.e. photons with energies larger than 13.6 eV), there will be a low density 
layer of ionized hydrogen at the surface of the disks. In the following, we ignore this layer as 
its contribution to the disk mass is negligible. Underneath this H II layer, there will be a layer 
dominated by far ultraviolet photons from the star where hydrogen is mostly neutral. As the 
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optical depth becomes large enough to shield this radiation, molecular hydrogen dominates. 
Even though molecular hydrogen dominates the disk mass in early phases of disk evolution, 
it is difficult to detect; due to the lack of dipole transitions, its ro- vibrational lines are very 
weak and the line-to-continuum ratio in the near-IR limits the sensitivity. Atomic hydrogen 
has electronic transitions, which need typical excitation temperatures in excess of 5000 K 
and therefore originate only in the surface layers very close to the central star, where the gas 
can be that hot. In this paper we study the role of the H I 21 cm line — a hyperfine transition 
of atomic hydrogen — for detecting the gas in various evolutionary stages of protoplanetary 
disks. 

The H I line could be an important tool to study the surface layers of young protoplan- 
etary disks. By simultaneous observations of the Hi 21 cm line and ro-vibrational lines of 
molecular hydrogen, we can observationally determine the H1/H2 fraction and compare it to 
predictions from disk models. This constrains both the photochemistry in the disk surface 
and the strength of the dissociating stellar (and surrounding) radiation field. The H I line is 
also a very important tool in traci ng disk dispersal and in particular photoevaporation pro- 



cesses. 



Hollenbach fc Gorti I (120061 ) have proposed that the FUV photoevaporation from the 
star itself is one of the main processes of disk dispersal. The FUV radiation can penetrate 
much deeper in to the disk than the EUV radiation and it heats the gas up to temperatures 



of a few 100 K. lAdams et al. I (120041 ) have shown that photoevaporation can already start at 

temperatures much lower than the virial temperature, T g as ~ 0.2 0.5 T^j-jt. The critical 

temperature is given by 

T clit = GM*m u /(kr) K, (1) 

where G is the gravitational constant, M* the stellar mass, mn the mass of an hydrogen 
atom, k the Boltzmann constant and r the radial distance from the star. Typical values 
for T crit are in the range of 180 to 10 000 K for our T Tauri disk models. This means, that 
the FUV radiation field can drive a photoevaporative wind from layers much deeper in the 
disk than the EUV and hence this wind would be more massive than the EUV wind and 
dominate the mass loss of the disk. Hi observations could thus be used to trace this mass 
loss and to derive photoevaporation rates, which can be compared to the theory. 

The rest of this paper is organised as follows. The Australia Telescope National Facility 
(ATNF) and the Giant Meterwave Radio Telescope (GMRT) H I 21 cm line observations are 
described in Sect. |2j In Sect. [3] we derive approximate estimates of the fraction of atomic 
hydrogen in circumstellar disks at different stages of evolution. Detailed models of each 
evolutionary stage ranging from protoplanetary disks to debris disks are presented in Sect. @] 
while in Sect. |5] we outline the details of the Hi 21 cm line radiative transfer. In Sect [6] we 
compare the model predictions with the observations, and finally in Sect. [7] we assess the 
prospects of using future radio telescopes to observe H I 21 cm line emission to trace gas in 
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circumstellar disks at various evolutionary stages. 



2. Observations 

Atomic hydrogen 21 cm line emission has never been detected in circumstellar disks. 
The only case for which upp er limits are avai l able in the literature is (3 Pictoris, which has 



been observed at ATCA by iFreudling et al. I (119951 . F95). This program was set up as a 



pilot survey. Our sample was hence selected to include disks from which HI was likely to 
be detectable (based on the presence of detectable amounts of molecular gas), while at the 
same time trying to cover a wide range of stellar and disk types. 

We have carried out 21 cm line observations of the disks around HD 135344, LkCal5 
and HD 163296 at the Giant Metrewave Radio Telescope operated by the National Centre for 
Radio Astrophysics, India. The stellar and disk classification and properties are summarized 
in Tabled! The observations were conducted on August 29 and 30, 2004. The total integra- 
tion time on each source was 4 hours. We used a 1 MHz bandpass with 128 channels, which 
results in a channel width of 1.6 km/sec. The standard calibrators 3C 232 and 3C 286 were 
observed at the start and end of the observing run and used to calibrate the visibility am- 
plitudes and the bandpass shape. Phase calibration was carried out with continuum sources 
1522-275, 0431+206 and 1751-253 for HD 135344, LkCal5 and HD 163296, respectively. For 
each disk, separate maps were produced at angular resolutions of 2"x3" and 7"x9" using 
appropriate UV ranges and tapers. The spectra at the position of the disks were in each 
case closely inspected for any signs of H I emission. Finally, all spectra were also smoothed 
with a 3 channel boxcar function. No evidence of any line emission was found in either the 
smoothed or unsmoothed spectra. One a upper limits are listed in Table [2J 
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Table 1: Observed stellar and disk parameters: Mdi S k(l-3mm) and Mdi S k(CO) are the total 
disk masses as derived from 1.3 mm continuum and CO radio observations respectively. 



star 



classification 



d 

[PC] 



T e ff 
[K] 



R* 
[R Q ] 



[M ] 



M, 



disk(1.3mm) 



M, 



disk(CO) 



1(T 2 M 10~ 2 M [erg s" 



j3 Pictoris A5V debris disk 

HD 163296 AlVe pre-MS 

HD135344 F4Ve pre-MS 

LkCal5 K5 pre-MS 



19.3 1 
122 2 
84 2 
140 4 



6 
65 3 
60 3 
60 5 



8200 1 
9230 2 
6200 2 
3980 4 



1.47 1 

2.7 2 

1.3 2 



1.75 1 0.003 4 
2.4 3 6.5 4 
1.3 3 0.28 4 
3.3 4 



2.6 x 10 25 7 
0.056 4 4 x 10 29 8 
0.00021 4 ... 
0.014 4 < 4 x 10 29 9 
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Table 2: Observed 21 cm line emission upper limits 



star 


beamsize 


Av 


la limit 


source 




[" x "] 


[km/sec] 


[mJy/beam] 




j3 Pictoris 


27x 25 


1.7 


3.7 


F95 




27x 25 


8.3 


1.7 


F95 




68 x 65 


1.7 


10.8 


F95 




68 x 65 


8.3 


4.7 


F95 


HD163296 


3x2 


1.6 


3.2 


this work 




9x7 


1.6 


13.7 


this work 




3x2 


4.8 


2.0 


this work 




9x7 


4.8 


10.1 


this work 


HD135344 


3x2 


1.6 


4.6 


this work 




9x7 


1.6 


12.5 


this work 




3x2 


4.8 


2.6 


this work 




9x7 


4.8 


10.6 


this work 


LkCal5 


3x2 


1.6 


7.6 


this work 




9x7 


1.6 


20.0 


this work 




3x2 


4.8 


5.1 


this work 




9x7 


4.8 


11.6 


this work 
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3. Estimates of the Hi mass fraction 

In the following we do not attempt to generate disk models that fit all the observed 
stellar and disk properties (see Tabled]), but rather use generic representative disk models 
for T Tauri disks, Herbig Ae disks, transition phase and debris disks to interpret the observed 
upper limits and to help plan future observations. 



Detailed chemo-phy sical disk models (IKamp fc Dullemond Il2004r iJonkheid et al. 112004 



Nomura fc Millar 1120051 ) have shown that the chemical structure of the surfaces of UV dom- 
inated protoplanetary disks resembles that of photon dominated regions seen at the surfaces 
of molecular clouds. The surface layers are dominated by far ultraviolet photons (6 - 13.6 eV) 
and contain mostly atomic hydrogen, and other atomic and ionized species. The transition 
to molecular species and thus to H 2 occurs roughly around a UV continuum optical depth 
r(1000 A) ~ 1. This transition occurs very close to the disk surface and thus the mass in 
this surface layer will naturally be small compared to the total disk mass. At very high gas 
temperatures (few 1000 K) however, the H/H 2 balance is dominated by H 2 destruction via 
collisions with neutral oxygen. 



In the case of you ng; active stars with high X-ray luminosity (~ 10 28 ' 5 — 10 31,5 erg s 1 



Glassgold et al. 1 119971 ). the chemistry and thus the depth of the surface layer is dominated 
by X-rays. Close to the star, X-ray heating dominates the disk energy balance, thus driving 
the chemistry that sets the H/H2 balance. In the outer disk, secondary ionization from X- 
rays determines the optical depth at which the transition from atomic to molecular hydrogen 
occurs. 

Before we go into the details of the complicated chemo-physical disk models in the next 
section, we outline the basic physics governing the H/H2 balance and hence the Hi mass 
fraction in the various evolutionary stages from young protoplanetary to very old debris 
disks. 



3.1. Flaring protoplanetary disks 

Where X-rays can be neglected, the hydrogen chemistry is driven by UV irradiation. 
We can estimate the H I mass in the disk by assuming that the optically thin surface layer of 
the disk contains only atomic hydrogen and that the transition to molecular hydrogen occurs 
around Ti n (1000 A) ~ 1. This might slightly overestimate the mass of atomic hydrogen as 
we do not take into account H 2 self-shielding. However, we do take into account the flaring 
angle of the incoming radiation. Given the UV absorption cross section of our disk models, 
cr(1000 A) = cr(UV) = 5.86 x 10 -22 cm" 2 (H-atom)" 1 and a typical flaring angle of a = 0.05 
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(appropriate for the typical luminosity of T Tauri stars, L* = 1.4 L Q ), we obtain a total 
column density of 

rv 

N tot = - = 8.5 x 10 19 cm" 2 (2) 
a 

for the optically thin surface layer. Assuming that the disk extends from 5 to 300 AU, this 
translates into a total mass of 2 x 10~ 6 M . Adding mass to an optically thick disk simply 
adds mass to the interior, while the surface layer remains mostly unchanged. This makes 
the H/H 2 fraction of the disk diminish with increasing disk mass (Fig. [1]). Herbig Ae stars 
are generally more luminous and hence their typical flaring angle is up to 50 % larger than 
that of the T Tauri stars. This increases the H/H2 mass fraction by the same amount. 

The pure UV case presents a lower limit to the estimated Hi mass in the disk. In 
the presence of X-rays, the column d ensity of atomic hydrogen in the disk surface increases 



substantially. iGlassgold et al. I (120061 ) assumed a typical X-ray luminosity for T Tauri stars 
of 2 x 10 30 erg s _1 . They find a total column density of iV tot ~ 8.5 x 10 21 cm -2 for the 
Hi surface. This number is fairly constant over the 5-20 AU range and hence we use it as 
an upper limit for the H I mass increase in the presence of X-rays. Such X-ray luminosities 
may lead to a significantly larger H/H2 mass fractions compared to the pure UV case. The 
shaded area in Fig. [T] indicates thus the possible range of atomic hydrogen masses ranging 
from the pure UV to the UV plus X-ray case. 




Fig. 1. — (a) Mass fraction of neutral hydrogen as a function of total disk mass. Flaring 
disks are shown as solid lines, whereas dashed lines are used for debris disks. The filled circle 
denotes the H I mass fraction in the HD 141569A disk model and the open circles denote the 
respective values for the series of (3 Pictoris models. Note that the (3 Pictoris models all have 
the same dust mass, the gas-to-dust mass ratio therefore differs for different models. The 
shaded area illustrates the range of atomic hydrogen fraction in the presence of X-rays, (b) 
Same for total H I mass. 
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3.2. Transition phase disks 



Massive protoplanetary disks are believed to evolve through a transition phase into 
debris disks. These transition phase disks are usually characterized by an absence of near- 
IR excess and thus large inner holes. The outer disks however, are still optically thick. 
Observations of the very few known transition pha se disks reveal that they generally have a 



very complex structur e (e.g. rings, spirals, clumps iMouillet et al. Il200ll ; iGrady et al. 
Fukagawa et al. 112004 ). which is difficult to treat analytically. 



2001 



We will use later i n this paper a disk model for the transition phase object HD 141569A 
( iJonkheid et al. 1120061 ) to fill the gap between the optically thick young disks and the opti- 
cally thin debris disks. This model has been constructed to fit the dust scattered light and 
CO line observations of HD 141569A. 

An UV only model predicts a relatively large Hi flux. In fact, Fig. [T] shows that the Hi 
mass fraction peaks around the regime of transitional objects such as HD 141569A. This is 
partly due to the fact that such transition phase disks are more tenuous than massive disks 
and hence UV photons penetrate deeper. 



3.3. Debris disks 

In very low mass disks, the H I mass fraction becomes 100% as the disk becomes totally 
optically thin even to UV line photons. However, such disks are very unrealistic, because 
generally dust grain growth and gas dispersal go along with a decreasing disk mass. Hence, 
there is a natural transition from protoplanetary disks to debris disks. 

In the next section, we compute debris disk models which apply to the UV only case 
and the predicted H I presents a lower limit for very active stars such as AU Mic, which have 
a significant X-ray luminosity. These models use gas-to-dust mass ratios varying between 
100 (i.e. value typical of molecular clouds) and 2. This affects the formation of molecular 
hydrogen as the increasing dust surface area per gas changes the gas thermal balance and 
hence the temperature dependent H 2 formation rate. The H 2 formation time scale is generally 
faster than 1000 yr even in the very low mass debris disks and equilibrium chemistry is thus 
a valid assumption. 

In debris disks, we generally have a size distribution of dust grains with a minimum 
and maximum grain size, a min and a max . However, the H 2 formation depends to first order 
only on the grain surface per hydrogen atom n^aa. A distribution of grain sizes, can hence 
to first order be replaced by grains with a fixed size a, where a is given by a 2 = a m i n a max 



( IKriigel 



2003 



, p. 276). In the following, we have chosen a fixed grain size a = 3 /xm. 



Since the dust surface area per hydrogen atom n^a^ stays constant in our models, the 
main impact of a decreasing gas-to-dust mass ratio is an increasing gas temperature and 
thus larger thermal velocities v. This increases the collision rate between hydrogen atoms 
and grains and hence the H2 formation rate. The photodissociation of H2 is not affected by 
the gas-to-dust mass ratio. So, the H/H 2 balance shifts towards molecular hydrogen as the 
gas-to-dust ratio decreases. We use the detailed results from our models described in Sect. @] 
to approximate the Hi mass fraction as a function of the gas-to-dust mass ratio 5. Fig. [T] 
shows that the Hi mass fraction scales approximately as 5 ' 7 (dashed line) for the debris 
disks. 



4. Detailed Disk Models 



This section describes the numerical models which we construct for the circumstellar 
disks. In a later section we use these models to predict the Hi 21 cm line emission. We 
use two different types of disk models here for the optically thin and thick systems. Our 
models do not include X-rays and deal with the pure UV case. The debris disks are typically 
optically thin and thus we solve the chemistry and heating/ cooling along radial rays starting 
in the top disk layers and moving down to the midplane. In the case of optically thick models, 
we divide the disk in vertical slices and calculate the coupled chemistry and heating/cooling 
in each slice. The different geometries are shown schematically in Fig. [2j These two different 
approaches are necessary, because in the debris disk systems, the main source of energy 
input is the direct stellar irradiation, while in the younger optically thick flaring systems, 
the irradiation is scattered from the surface into the disk. In the following, we summarize 
these two disk modeling approaches and the free parameters entering the disk models. 
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Fig. 2. — Geometries and illumination for the (a) debris disks and (b) flaring disks. The 
dotted lines indicate how the models discussed in Sect. H]are split up in slices. The debris 
disks are calculated from top to midplane following radial cuts, the flaring disks from top to 
midplane following vertical cuts. 
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4.1. Disk physics and chemistry of optically thin disks 



Th e basic physics and chemistry of these models is described in iKamp fe Bertoldi 



( I2000l ). IKamp fe van Zadelhoff I (120011 ). Here we only briefly summarize the salient features 
of these models. The disk density structure is assumed to be, n(r,z) m r _p exp — z 2 /2h 2 , 
where we further assume that the scaleheight h scales linearly with radius, and h/r = 0.15. 
The radial power law exponent p is set to 2.5. The chemistry is modelled using a network 
consisting of 48 different species covering the elements H, He, C, O, S, Mg, Si, and Fe. We 
use typical IS abundances for these elements (Table Ej). These species are connected through 
281 reactions, including also cosmic ray chemistry, photochemistry and the chemistry of 
excited H 2 . The dust temperature is derived from the assumption of large black body grains 
in radiative equilibrium around the star, a = 3 /im, p = 3.0 g cm -3 . The gas temperature 
is derived from a detailed energy balan ce including the most relevant heating and cooling 
processes (IKamp fe van Zadelhoff 11200 ll ). Table 0] lists the basic parameters of the computed 
disk models. 
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Table 3: Element abundances in the disk models normalized to a hydrogen density of 
10 12 cm~ 3 for debris disk models 



element 


abundance 


element 


abundance 


H 


12.00 


S 


6.90 


He 


11.00 


Mg 


6.00 


C 


8.15 


Si 


5.90 





8.51 


Fe 


5.40 
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H2 photodissociation occurs in the wavelength range between 912 and 1110 A. We 
calculate the photodissociation rate from 



T = 4.2 X 10- U (x/shield,r 



G f s hieid,v)n(H 2 ) s 1 cm 3 



(3) 



where x is t he integrated s tellar radiation field between 912 and 1110 A nor malized to the 
Habing field (jHabing Ill968l ). The stellar radiation field is an ATLAS9 model (IKurucz Ill992l ) 
with stellar effective temperatures as listed in Tab. HI An interstellar radiation field with 
Go = 1 in units of the Habin g field is added to t he st ellar one. The H 2 self shielding 
is approximated by Eq.(37) of iDraine fc Bertoldi I (119961 ) for radial (/ s hieid,r) and vertical 
directions (/ s hieid,v)- We assume here that radiative decay of vibrationally excited H 2 is 
faster th an direct ph otodissociation out of these levels. This approximation is valid for 
X < 10 5 (jShull Ill978l ). a limit that holds for all models, except the inner 5-10 AU of the T 
Tauri disk and the inner 1.4-2 AU of the Herbig Ae disk. Even though the H/H 2 fraction 
in those areas might change if direct photodissociation from vibrationally excited levels is 
included, the mass of those inner disks is negligible compared to the total disk mass. The 
predicted H I line profiles and fluxes will hence not be significantly affected. 



The H 2 formation rate is 



R 



form 



0.5nnnd<JdVe-R 2 S(T) s cm 



.-3 



(4) 



where tt-h is the hydrogen number density, n& is the dust grain number density, o"d is the dust 
grain cross section, v the thermal velocity, en 2 the recombination efficiency of hydrogen atoms 
on the surface and S(T) a temperature dependent sticking coefficient for atomic hydrogen on 
the grain surface. Higher gas temperatures lead to larger thermal velocities and hence to more 
H-grain colli sions. The recombination efficiency depends on the details of the H 2 formation 
mechanism. ICazaux fc Tielens I (120021 ) approximate the recombination efficiency 6j^2 as a 
function of grain temperature. While at low grain temperatures nearly every incoming 
H-atom recombines with a surface H-atom to form H 2 , this process becomes much less 
efficient at higher temperatures (50-100 K), where direct re-evaporation of pysisorbed H- 
atoms competes with H 2 formation. The sticking probability S(T) can in principle depend 
on temperature, but here we assume it to be equal to one. 

The models for (3 Pictoris are calculated assuming a cons tant dust mass of 1 .3 10~ 6 M 
derived from fitting the infrared spectral energy distribution (jChini et al. Ill99ll ). Since the 
total gas mass in the disk around this star is not yet well constrained, we assume here different 
gas-to-dust mass ratios ranging from the canonical value of 100 down to 2. Fig. [3] shows the 
H and H 2 abundances for two vertical cuts at 50 and 200 AU in the three disk models. This 
quantitatively illustrates what has been stated in Sect. 13.31 the gas temperature increases 
as the gas-to-dust mass ratio decreases. The heating/cooling balance shifts, because the 
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line cooling decreases with decreasing gas densities, while the photoelectric heating stays 
constant. 



r=50AU n . r=200AU ,. 

'I -o I U 

5 10 15 20 25 20 40 60 80 100 

z [AU] z [AU] 




5 10 15 20 25 20 40 60 80 100 

z [AU] z [AU] 



Fig. 3. — H/H 2 transition for two vertical cuts at 50 and 200 AU in the three debris disk 
models: (a), (b) 44 M Earth , (c), (d) 8.8 M Ea rth, (e), (f) 0.88 M Earth . The solid lines are the H 
abundance, the dashed lines the H 2 abundance. The dotted lines denote the gas temperature 
using the scale on the right axis. 
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Table 4: Parameters of disk models: debris disk models, HD 141569A is a transition phase 
disk model and flaring disk models for a Herbig Ae and T Tauri star 



type 


Toff 


R* 


M* 


Mdisk 




M gas /M dust 


Rin 


Rout 


a 


cr(UV) 


Luv 




[K] 


[R©] 


[M ] 


[M ] 






[AU] 


[AU] 


H 


[cm 2 /(H-atom)] 


[erg s" 1 ] 


debris disk 


8250 


1.7 


2.0 


1.3 x 10" 


-4 


100 


40 


500 


3 


2.34 x 10~ 23 


1.9 x 10 27 










2.6 x 10" 


-5 


20 


40 


500 


3 


1.17 x 10~ 22 


1.9 x 10 27 










2.6 x 10" 


-6 


2 


40 


500 


3 


1.17 x 10~ 21 


1.9 x 10 27 


HD 141569A 


10 000 


1.7 


2.0 


2.4 x 10" 


-4 


100 


80 


500 


0.1 


5.86 x 10~ 22 


2.7 x 10 33 


Herbig Ac 


9500 


2.4 


2.5 


0.01 




100 


1.4 


300 


0.1 


5.86 x 10~ 22 


6.6 x 10 29 


T Tauri 


4000 


2.5 


0.5 


0.01 




100 


5.0 


300 


0.1 


5.86 x 10~ 22 


2.0 x 10 31 
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4.2. Disk physics and chemistry of optically thick disks 



For the massive optically thick disk models we follow a similar approach as for the 
optically thin models. However, the underlying equilibrium density structure is very different ; 



it is taken from two-dimensional hydrostatic equilibrium models (IDullemond et al. II2002I ). 



In these models, the gas temperature is first approximated by the dust temperature found 
from detailed continuum radiative transfer. We then re-evaluat e the gas temperature in th e 
context of the coupled chemistry and heating/cooling balance (IKamp fc Dullemond 112004 ). 
Another difference from the hydrostatic equilibrium models is need to account for the flaring 
angle, that is the grazing angle under which the stellar irradiation hits the disk surface (see 
Fig. ED. 

The stellar radiation field for a T Tauri star is a combination of an ATLAS9 model and 



a scaled solar chromosphere model as outlined in IKamp &: Dullemond I (12004) . Even though 
the chromospheric UV excess contributes little to the total luminosity and hence the flaring 
of the disk, it is relevant for the photodissociation of molecules and the photoelectric heating 
of the disk. In the case of the Herbig Ae star, we used a pure ATLAS9 photospheric model 
with T e g = 9500 K and logg = 4.0. Fig. 0] compares the various UV and optical radiation 
fields used in the disk models. Diffuse interstellar UV irradiation is not relevant for young 
disks at radii smaller than 300 AU. 

The chemical network and heating/ cooling processes are the same as for optically thin 
disks. The dust properties however are different reflecting a much younger evolutionary stage. 
We use the same gr ain properties and opaciti es as used in the underlying continuum radiative 



transfer models by IDullemond et al. I (120021 ) , a = 0.1 /mi, cr(UV) = 5.86 x 10 22 cm 2 /H- 
at om The protoplaneta r y disk mo dels described in Tab le H] have been previously published 



in 



Kamp fc Dullemond I (12004 ) and IKamp et al. I (120061 ) as a typical T Tauri and Herbig Ae 



star respectively. 



5. Hi 21 cm line radiative transfer 

We assume here that the level population numbers are in LTE, i.e. T ex = T gas , which is 
a valid approximation for nn 4 10~ 4 T gas . Given typical particle densities and gas temper- 
atures of 10 3 < nn < 10 11 cm -3 and 10-10 000 K, this assumption is fullfilled throughout the 
entire disk. We thus assume a ratio of rii/n = 3 for the level population numbers, a good 
approximation for T gas > 10 K (error smaller than 1%). Even though the gas temperature 
does not play a role in the Hi level populations, it is still crucial for the determination of 
the H/H2 chemical balance and hence the total number density of atomic hydrogen. 
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Fig. 4. — UV radiation fields for the various disk models: ATLAS9 (3 Pictoris model (solid 
line), ATLAS9 Herbig Ae star model (dotted line), T Tauri model with UV excess (dashed 
line) . 



The line emission and absorption coefficients eio and Kq\ are derived from 

e 10 = fi,i/ 10 Alotii(H) erg s" 1 cm" 3 (5) 

«oi = — ^5oi™o( H ) s cm (6) 

The absorption coefficient is corrected for stimulated emission assuming LTE level population 
numbers and using the first order term of the expansion of the exponential. The basic atomic 
parameters for the line are A 10 = 2.87 fO" 15 s" 1 , v w = f.4205 GHz, and ni(H) ~ 3/4 x n(H). 
The total emission per velocity bin dv is derived from an integration along the optical path 
within the beam area dA and within the respective velocity bin dv =[vi,Vi + i] 



'10 



/ — ^- exp(— drio) ds erg s 1 km s 1 . (7) 

J beam 



The optical depth drw is derived from the absorption coefficient 

dr w = I ds . 



dv 

The typical 'beam averaged' optical depth are smaller than r ~ O.OOf for the Herbig Ae 
disk seen edge-on and ~ 0.02 for the pole-on model. Even though the column densities are 
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smaller in the pole-on model, the velocities in the line of sight are zero; hence the optical 
depth contribution at v — is maximized. 

The optical depth is - in the case of a rotating disk - defined as a function of velocity, 
t(v). Hence, the corresponding column density to derive the optical depth is the averaged 
column density of all H-atoms falling in the velocity bin dv averaged over the beam. We 
assume that the disk is in keplerian rotation v^ ep around the star with central mass M*. The 
velocity profile for an inclined disk (the inclination angle % is defined between the line of sight 
and the disk midplane, so that i = denotes an edge-on disk ) is given by 

v TOt (r, z, 0) = v kep (r, z) sin(0) cos(?) , (9) 

where is the azimutal angle in the disk. 

To calculate observed 21 cm line fluxes, the disk is projected onto a Cartesian grid with 
the z-axis in the direction of the line of sight. A separate cube is created for each velocity 
bin (see Fig. 5). For each of the cubes, the emission and optical depth are integrated along 
the line of sight. 
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Fig. 5. — Velocity binning (a) and coordinate transformation (b) for the 21 cm line radiative 
transfer. The shaded area in (a) corresponds to the velocity bin [v — dv/2,v + dv/2] and 
matter within that area will contribute to the line emission profile at the velocity v. 



The intensity arriving at Earth is given by 



I 



10 — 



4nd 2 



1 



/ 



F w dA erg cm 



2 s 1 km s 



-i 



(10) 



Converting to more common units of erg cm 2 s 



1 Hz 1 is done by using 



dv 

dv = c — . 
Via 



(11) 



We use the beamwidth given in Table [2] as a FWHM to convolve the resulting images in each 
velocity channel assuming a gaussian beam profile g(x,y). This leads to the final intensity 
in observers units, mJy per beam, extracted at the center (xi, y\) 



We will constrain disk parameters, especially the total disk gas mass, by comparing the 
observed upper limits for the Hi 21 cm line presented in Sect. [2] to predicted emission from 
the individual disk models. In Sect. [3] we presented rough estimates of the Hi mass in the 
various disks based on various simplifications. In this section we take a more quantitative 
approach involving detailed disk models as well as stellar parameters, inclinations and dis- 
tances appropriate for the individual sources. We apply the models described in Sect. H]to 
predict H I line maps and intensities based on the observed beam sizes. Table [5] summarizes 
the basic H I characteristics of those disk models. 




mJy/beam . 



(12) 



6. 



Predictions for H I disk emission 
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Table 5: Hi 21 cm line characteristics of the disk models 



type Mdisk H I mass beam size peak intensity peak surface brightness 

[M ] [M Earth ] [AU] [mJy] [^Jy/arcsec 2 ] 



debris disk 


1.3 x 10" 


-4 


4.0 


1280 


18.0 


480 




2.6 x 10" 


-5 


0.4 


1280 


1.8 


40 




2.6 x 10" 


-6 


9 x 10~ 3 


1280 


4.4 x 10- 2 


1 


HD 141569A 


2.4 x 10" 


-4 


13.0 


1100 


1.6 


290 


Herbig Ae 


0.01 




0.3 


980 


0.05 


110 


T Tauri 


0.01 




1.3 


1120 


0.24 


400 
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6.1. (3 Pictoris 

The disk around (3 Pictoris is the best studied examp le of a debris disk. O VI line 



emission has been detected with FUSE (IDeleuil et al. Il200ll ) indicating some level of stellar 



activity. A followup search for the O VII line wit h XMM- Newton has revealed a very low flux 



of 6 x 10~ 16 erg cm" 2 s -1 f lHempel et al. II2005I ). This results in a X-ray to UV luminosity 
ratio of 10~ 2 . The impact of such an X-ray flux on the chemistry of debris disks (gas-to-dust 
mass ratios smaller than 1) is yet unknown. We can therefore not rule out the possibility 
that the disk around (3 Pictoris contains additional atomic hydrogen produced by X-rays. 

The dust mass in the /3 Pictoris disk is well known from previous observations and we 
assume in this paper a fixed dust mass of 0.44 MEarth (see Sect. 14. ip . We computed models 
with gas-to-dust mass ratios from 100 to 2. Fig. [6] shows the distribution of atomic hydrogen 
in the disk for individual velocity channels (Av = 1.0 km/s). The atomic hydrogen is located 
in a thin layer on top of the entire disk. The disk towards the midplane looks filled by atomic 
hydrogen, but this is due to the projection of the three dimensional disk model into the plane 
of observations. It stems from H I in the disk surface outside of the observing plane. This 
can be seen in Fig. [7J where we plot the integrated Hi 21 cm line image on top of the 
WFPC2 dust scattererd light image of Al Schultz and his team (private communication). As 
expected, the maximum of the Hi emission occurs at the surface of the gaseous disk. Note 
also that the HI emission lies well within the region from which scattered light from dust 
grains can be detected. Recent imag ing of the disk in the Nal D lines indicates scaleheights 



of h/r ~ 0.3 (IBrandeker et al. 1120041 ). a factor 2 higher than what our models actually use. 
It is currently unclear which processes lead to such large scaleheights, but one reason might 
be higher gas temperatures in case of a very low gas-to-dust mass ratio. As explained in 
Sect. 14. H this is due to a shift in the thermal balance between photoelectric effect and fine 
structure line cooling. The distribution of the micron sized dust grains responsible for the 
scattered light is sustained by collisions among dust grains and/or larger bodies within the 
disk and is not coupled to the gas. 

We checked the resulting line intensities for the different beam sizes used in the obser- 
vations. The beams were always centered on the stellar position. It turns out that the 66.5" 
(1280 AU) beam includes most of the Hi emission of (3 Pictoris, which is not surprising given 
our disk model extension of r = 500 AU (see Fig. [BD- The smaller beams miss a large fraction 
of the emission at intermediate velocities (v ~ 2 km/s), which is located in the outer regions 
of the disk (Fig. [6]). Given the density law of the underlying disk model, most of the disk 
mass — and hence also the H I mass — resides in the outer parts of the disk. The maximum 
H I intensity of the 44 MEarth disk model, 18 mJy, is larger than the observed upper limit of 
4.7 mJy. The ATCA observations (F95) therefore imply an upper limit of ~ 15 MEarth on 
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the total gas mass (H1+H2) in the disk around (3 Pictoris (Fig. [9]). 
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Fig. 6. — Map of the atomic hydrogen distribution in the (3 Pictoris 44 M Earth disk model 
(gas-to-dust mass ratio of 100, inclination 0°). The different gray shades denote steps of a 
factor 10 in emission. 
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Fig. 7. — Hi 21 cm line image from the 4.4 Msarth mass model (blue) on top of the HST 
optical scattered light image taken with Wide Field Planetary Camera 2 (red, credit Al 
Schultz and his team). The maximum Hi emission (white) comes from the surface of the 
gaseous disk; peak fluxes are 17/zJy/arcsec 2 , well within the reach of the SKA. 
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Fig. 8. — 21 cm line intensity for the f3 Pictoris 44 M Earth disk model (gas-to-dust mass 
ratio of 100, inclination 0°) and four different beam sizes: 72.5" (1400 AU, solid line), 57" 
(1100 AU, dotted line), 33.7" (650 AU, dashed line), and 22.3" (430 AU, dash-dotted line). 
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Fig. 9. — 21 cm line peak intensities for different gas-to-dust ratios in the f3 Pictoris model: 
100 (M gas = 44 M Earth ), 20 (M gas = 8.8 M Earth ), and 2 (M gas = 0.88 M Earth ). Line fluxes were 
calculated with a gaussian beam with FWHM of 66.5" (1280 AU). The dashed line denotes 
the upper limit derived from the ATCA observations (Table [2]) for this beamsize. 
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6.2. LkCal5 



For modeling th e Hi e mission from this star, we used the T Tauri disk model by 
Kamp &: Dullemond I (120041 ) . This model exte nds from 5 to 3 00 AU and has a mass of 
0.01 M . The distance to the star is 140 pc and lQi et al. I (120031 ) have shown from modeling 
CO and HCO + that the inclination is 57° ±5°. LkCa 15 has not been detected in the ROSAT 



All Sky Survey (RASS). iNeuhauser et al. I (119951 ) put an upper limit of 4 x 10 29 erg s" 1 on the 
X-ray luminosity for this source. The resulting upper limit to the X-ray to UV luminosity 
ratio is 0.02. 

The Hi observations with a beam of 8" comprised the entire disk as opposed to the 
smaller beam of 2.5 , which contains only the inner 125 AU of the disk. Using the larger 
beam, we predict a 21 cm line flux of 0.24 mJy, two orders of magnitude lower than the 
upper limit derived from GMRT observations. For the smaller beam size, we obtain fluxes 
of 0.13 mJy. The observed upper limit of 20 mJy to the 21cm flux corresponds to an upper 
limit of 3 x 10~ 4 M for the atomic hydrogen mass. 

Since we have not accounted for the Hi production by X-rays, in young active stars 
which are strong X-ray emitters, the Hi mass would be much higher than what we obtain. 
As can be seen in Fig. [TJfor such stars in nearby (50-100 pc) young associations, the 21 cm 
line flux might be detectable. 



6.3. HD 163296 



For mod e ling t he Hi emission from this star, we used the Herbig Ae disk model by 



Kamp et al. I (120061 ). The stellar effective temperature in this model of 9750 K is only 
marginally higher than the value quoted in Table H] an d the stellar mass is 2 .5 M Q . HD 163296 
is a weak X-ray source with L x = 4 x 10 29 erg s -1 (ISwartz et al. 1120051 ). corresponding to 
an X-ray to UV luminosity ratio of 0.6. Therefore, the predicted Hi mass from our model 
should be considered as a lower limit. 

The model extends from 1.7 to 300 AU and has a mass of 0.01 M^. We use a distance of 



122 pc and an inclination angle of 65° from I Jayawardhana et al. I (120011 ) and lDominik et al. 



( 120031 ) respectively. As in the case of LkCal5, we used the larger beam of 8" , so that the 
beam contains the entire disk. The peak Hi 21 cm flux for this model is 0.05 mJy. The 
fluxes are lower than for LkCal5, because the Hi fraction in the model is about an order of 
magnitude lower (~ 0.01 %). This is due to a weaker UV radiation field. While we assumed 
a chromospheric UV excess for the T Tauri star, we chose a pure photospheric model in the 
case of the Herbig Ae star; chromospheric activity is mostly restricted to solar-type stars. 
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The observed upper limit of 13.7 mJy corresponds to an H I mass upper limit of 3 x 10 4 M G 
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6.4. HD 141569A 

Our estimates of atomic hydrogen fractions in disks suggested that the fraction would 
be maximal in transiti on phase objects. HD 141569A which has been modeled extensively by 
Jonkheid et al. I (120061 ) is such a transition phase disk. This star has not yet been observed 
in the 21 cm line, but given that it is an excellent candidate for future observations, we 
present here predictions for the Hi 21 cm flux. 

The param eters for the du s t phas e of the model have been derived from scattered light 
observations of iMouillet et al. I (120011 ). while the gas phase constraints come from a fit to 
the CO rotational emission (J=3-2) from lDent et al. I (120051 ). The dust is mostly confined in 
two rings located at 185 AU and 325 AU. The gas follows a smooth surface density profile 
£ r -!- 2 extending all the way in to 80 AU. This star has two M-type companions at 
distances of 790 and 670 AU that are most likely the source of the X-ra ys measured in the 



ROS AT All Sky Survey. The X-ray luminosity is 8.72 x 10 erg s ( 1 Weinberger et al. 



20001 ). Given the distance of the X-ray source from the disk around the primary star, the 
X-ray flux is a factor 5 x 10 5 lower than the UV flux at 40 AU. Thus, X-rays will not affect 
the chemical processes around HD 141569A. 

The maximum emission is 1.6 mJy using a beam diameter of 11" (1100 AU). With the 
current instrument sensitivity, this would translate into a 1 a signal. However, given the 
caveats of the above presented disk model, we cannot exclude that there might be more H I 
than anticipated from CO observations. A different gas density profile or an underestimation 
of the UV radiation field of this star may lead to an even higher fraction of atomic hydrogen 
and hence to a larger signal. Clearly, this is the most promising candidate to look for atomic 
hydrogen. 



Discussion 



In the following, we summarize the main results of the Hi modeling and discuss them 
in the light of present and future observing possibilities. 
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7.1. Gas in the disk around f3 Pictoris 



The gas mass upper limit from the Hi observations is ~ 15 MEarth- This amounts 



to a Hi midplane column density of 1.4 x 10 cm . However, the same model predicts 
H 2 , Oi, Ci, and Cll midplane column densities of 4.1 x 10 21 , 1.8 x 10 18 , 1.2 x 10 17 , and 
2.3 x 10 16 cm -2 . These values are much larger th an t hose derived from UV absorption 



measurements by iLecavelier des Etangs et al. I (120011 ) and iRoberge et al. I (120061 ) . The latter 
paper compares the elemental composition in the stable (3 Pictoris gas to the composition of 
the Sun, meteorites and comets; it leads to the conclusion that the gas in the disk around 
this star is extremely carbon overabundant, raising the interesting question of the origin of 
this peculiar chemical composition. 

Current Hi 21 cm line observations cannot rule out the presence of enough atomic 
hydrogen in the disk around (3 Pictoris to obtain a solar gas composition - in terms of 
hydrogen versus metals; metals here comprise all elements besides hydrogen and helium. 
More sensitive radio observations are needed to determine the true hydrogen content of this 
disk. Since current facilities such as the GMRT and ATCA are limited to sensitivities of the 
order of 1-2 mJy, we ha ve to wait for the Square Kilometer Array (SKA) to address this 
question (ISchilizzi II2004I ). The SKA will have a roughly 10 times larger collecting area than 
current telescopes and a spatial resolution of up to 1 mas for imaging. The sensitivity of the 
SKA is expected to be at least a factor of 100 above the one of current interferometers over 
a wide range of resolutions. 

Our models, however, suggest that the dominant form of hydrogen in debris disks is 
molecular and not atomic hydrogen (Sect. 13.31 and Sect. I4.ip . Thus the FUSE non-detection 
of H 2 i n the j3 Pictoris UV spectrum puts the most stringent upper limit on the disk hydrogen 
mass (ILecavelier des Etangs et al. 1 1200 lh . The H 2 column density from these observations 
is lower than 10 18 cm -2 ; this leads to an upper disk mass limit of 0.005 Msarth (0.4 Mmooii)? 
if we apply the same underlying disk structure we used in our debris disk models. 



7.2. Future Observations using Hi as a tracer 

Our models show that the non-detection of the Hi 21 cm line is not due to the absence 
of atomic hydrogen in these disks, but rather due to the weakness of the signal as compared 
to the sensitivities of current telescopes. However, our models also show that the SKA will 
be able to detect Hi 21 cm emission from protoplanetary disks around stars in nearby star 
forming regions. These observations would yield a number of interesting constraints on the 
properties and evolution of these disks. 
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In massive disks, the Hi 21 cm line traces the thin layer of the disk surface that is 
directly exposed to X-ray and UV irradiation; the H/H 2 mass fraction depends mainly on the 
strength of the radiation field impinging the disk surface. More stringent upper limits would 
hence give an independent check of the irradiation environment of the star. In addition, the 
velocity information in the line profile provides an independent measurement of the stellar 
mass, disk extension and evaporative flow from the disk surfa ce. Typical velocitie s for the 



evaporative flow are larger than 20% of the local sound speed (lAdams et al. 1120041 ). 

c s = VGM,/(0.5r) = 42 ( (—J km s' 1 . (13) 

Such an outflow should be detectable at the typical line resolutions of 1 km/s used in the 
radio regime. 

At later stages, e.g. in the transition phase disk around HD 141569A, the mass fraction 
of atomic hydrogen is of the order of 10%. This makes H I indeed a good probe of total disk 
mass in such transition phase disks. (See Fig. [lb for the dependence of H I mass on total 
disk mass). These disks are the best candidates for 21 cm line detection; our models indicate 
that existing telescopes may be able to detect emission from disks at distances smaller than 
100 pc. 

Debris disk systems often have a very low gas-to-dust mass ratio, which drives the 
hydrogen into its molecular form. In such debris disks, the H I fraction is very low, smaller 
than 10%. The integrated Hi flux is therefore significantly lower than for the transitional 
phase disk. This scenario therefore predicts hugely different total H I emission for the different 
types of disks which can in principle be tested observationally. 

Finally, we note that we have ignored Hi production by X-rays, which means that our 
models underestimate the Hi 21cm flux. Self consistent modeling of the X-ray chemistry 
as well as X-ray heating will hence probably lead to an even more optimistic prediction for 
the detectability of the 21cm emission from circumstellar disks by next generation radio 
telescopes. 
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constructive referee report that helped to improve the content and readability of the paper. 
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